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ABSTRACT 


The computational methods used to predict and optimize the thermal- 
structural behavior of aerospace vehicle structures are reviewed. In general, 
two classes of algorithms, implicit and explicit, are used in transient ther- 
mal analysis of structures • Each of these two methods has its own merits • 
Due to the different time scales of the mechanical and thermal responses , the 
selection of a time integration method can be a difficult yet critical factor 
in the efficient solution of such problems. 

Therefore mixed time integration methods for transient thermal analysis 
of structures are being developed. This proposed methodology would be readily 
adaptable to existing computer programs for structural thermal analysis. 

1. INTRODUCTION 

Over the last two decades , significant attention has been devoted to the 
development of lightweight, durable thermal protection systems (TPS) for 
future space transportation systems . Research programs are currently under 
way at the Langley Research Center to investigate various metallic TPS concepts 
[1]. One of the proposed candidates is the titanium multiwall tile (see [2] 
and references therein for a discussion). Early design procedures of the TPS 
concept involved both analytical and experimental studies. In particular, a 
degree of confidence has been established in the TPS concept due to the design 
studies by Jackson and Dixon [3] and Blair et. al. [4]. 

A titanium multiwall tile consists of alternating layers of superplas- 
tically formed dimpled sheets and flat septum sheets of titanium foil. As de- 
scribed in reference [3], this multiwall concept impedes all three modes of 

* 

The partial supppgrt of NASA under Grant No. NAG-1-210 to this research is 
gratefully acknoMedged. 


161 



heat transfer conduction, radiation and convection* The superplastlcally 

formed dimpled sheets and the long thin conduction path tend to minimize heat 
conduction. The flat septum sheets of titanium foil impede radiation* The 
small individual volumes created by the dimpled layers virtually eliminate air 
convection* The optimal design of such thermal protection systems requires 
effective techiques in coupled thermal and stress analyses* Finite element 
methods offer the greatest potential in modeling such complicated problems* 
However, the resulting semi -discrete equations may involve many thousand 
degrees of freedom* Since the problem to be solved is transient and non- 
linear, the selection of an appropriate time integration method is an essen- 
tial step in the solution of such a complicated problem* Adelman and Haftka 
[5] recently conducted a survey study on the performance of explicit and 
implicit algorithms for transient thermal analysis of structures* Calcula- 

tions were carried out using the SPAR finite element computer program [6] and 
a special purpose finite element program incorporating the GEARB and GEARIB 
algorithms. Based upon their studies, they concluded that, generally, implic- 
it algorithms are preferable to explicit algorithms for "stiff" problems , 
though non -convergence and/or wide-banding of the resulting matrix equations 
may decrease the advantage of the implicit methods* 

These difficulties are similar to those found in fluid-structure prob- 
lems. Over the past few years, several remedies have been proposed for these 
difficulties. Belytschko and Mullen [7J have proposed an explicit-implicit 
method where the mesh is partitioned into domains by nodes and the partitions 
are simultaneously integrated by explicit and implicit methods. Hughes and 
Liu [8J have proposed an alternate implicit-explicit finite element method 
where the mesh is partitioned into domains by elements and this element parti- 
tion concept simplifies the computer-implementation and enhances its compati- 
bility with the general purpose finite element software. 

Although the implicit-explicit method has been proven to be very success- 
ful in some fluid-structure interaction problems (see e.g. , [8-10] ) , the size 
and complexity of the program are increased because of the addition of the 
implicit method. To overcome these difficulties , Belytschko and Mullen [11] 
have proposed an E*^-E partition, in which explicit time integration is used 
throughout. However, different time steps within different parts of the mesh 
can be employed simultaneously. Partitioned and adaptive algorithms for ex- 
plicit time integration have also been proposed by Belytschko [12]. 

Recently, Liu and Belytschko [13] put forward a general mixed time 
implicit-explicit partition procedure within a linear context. It incorpo- 
rates the mentioned algorithms as special cases and is shown to have better 
stability properties than those in E^-E partition [11]. Similar concepts can 
also be used in transient conduction forced-convection analysis (see Liu and 
Lin [14]). 

In the present paper, we extend these implicit-explicit concepts^ nodes 
and elements^ to transient thermal analysis of structures where different time 
integration methods with different time steps can be used in each element 
group. The aim of this approach is to achieve the attributes of the various 
time integration methods . 
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For example, in transient structural analysis, explicit methods require 
the size of the time step to be proportional to the length of the shortest 
element, while in transient thermal analysis, explicit methods require the 
step size to be proportional to the the square of the length of the shortest 
element. So it is more advantageous to employ this mixed time implicit- 
explicit technique for transient thermal analysis of structures since the 

partition proposed in [11,12] is often inefficient for this kind of 
problem though it is very efficient in structural analysis. 

In section 2 we review the finite element formulation for transient heat 
conduction. In section 3 we describe the mixed time integration procedures 
viz two element groups ”A” and *'B” . A family of integration partitions can 
then be deduced by selecting the appropriate definitions for the quantities of 
**A” and "B“. Five useful partitions which are of practical importance are 
presented. In section 4 the stability characteristic of the algorithm is 
discussed. In section 5 we generalize the mixed time methods described in 
section 3 to NUMEG element groups. A computational algorithm for this mixed 
time implicit-explicit integration is also presented. Numerical results are 
presented in section 6 and conclusions and suggestions for further research 
are presented in section 7. 

2. FINITE ELEMENT FORMULATION FOR TRANSIENT HEAT CONDUCTION 

We consider a body ^ enclosed by surface T which consists of two parts: 
r and r • The Cartesian coordinates of the body will be denoted by x. . 

g q ^ 

The governing equations for transient heat conduction are : 



in 

(2.1) 

0 = g 

for X. in r 
i g 

(2.2) 


for X. in F 

i q 

(2.3) 

and 

o 

CD 

II 

CD 

for x^ in Q and t = 0 . 

(2.4) 

Here a comma designates a partial derivative with respect to x 

^ ; a super- 

script dot designates ^time (t) 

derivative; n. is the component of 
thermal diffusivity (the ratio of 

the outward 

unit normal vector; C is the 

thermal con- 


ductivity to specific heat times density); 0 is the temperature; h is the 
convective heat transfer coefficient; and g, q and 0^ are given functions. 
Repeated indices denote summations over the appropriate range. 

The variational or weak form of equations (2.1)-(2.4) is: 

(0,v) + a(0,v) == (q,v) in (2.5) 

where v is the test function; and 

(0,v) = / ^ 0 vdil (2.6) 
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(2.7) 


a(e,v) - e,^v,^dn + hev dr 


and 


(q,v)j, * /j, qv dr (2.8) 

q q 

The finite element equations are obtained by approximating the trial 
functions by shape functions (N^) so that 


V = 


and 


e 


NEQ 

lii 


(2.9) 

NUMNP 

Z 


(2.10) 

1=NEQ+1 

V + g 


(2.11) 


Here NUMNP is the total number of nodal points used in the finite element 
mesh and NEQ is the number of trial functions used (for this particular case 
it is equal to the number of equations to be solved). 


The resulting semidiscrete equation for transient heat conduction is 


then: 


Me + K 0 = F 


with initial condition 

e(0) = e 

' ~o 


( 2 . 12 ) 


(2.13) 


where 

M . [M^j] - Nj) - /a ^ Vjda (2.14) 

K - [K^j] - It “‘i"j >"■ «.15) 

and 

£ = [FiJ = (q.Njj. - a(Nj^,Njgj^ (2.16) 

q 

M and K aje assumed to be symmetric and positive definite. The thermal para- 
meters, C and h, are in general temperature dependent. However for discus- 
sion purposes herein, they are assumed constant throughout. 

3. MIXED TIME PARTITION PROCEDURES 

In this section, mixed time integration methods are employed to solve 
equations (2.12) and (2.13). For the purpose of describing these mixed time 
integration techniques we subdivide the mesh into element groups A and B, each 
of which is to be integrated by a different method. Let n be the time step 
number, 0 , v and F be approximations to 0(t ) , &(t ) and F(t ) re- 

spectlvely. Let mAt and At be the time steps used for element group A and 
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tlvely. Hence it follows that .any^gli 
trices, cf, M - and K - K^+ . 

/s/ 


element group B respectively, where m Is an Integer and is greater or equal to 
1. A time step cycle (mAt) can then be defined by an Increment of m substeps 
with a time step of At each, so that one time step cycle Is defined by step n 
to step n+m. The portions of the matrices obtained by assembling element 
group A and element group B are denoted by superscripts "A” and ”B”, respec- 

^global matrix Is the sum of the two ma- 
Nodes associated with only element 
group B are denoted by superscript ”B", whereas those which are in contact 
with at least one element of group A are denoted by superscript "A"; nodes 
which are connected to both group A and group B are designated by ”C", so ”C” 
is a subset of "A”. To simplify the presentation, we further denote those 

element^matrices assoclat^ with ^t least one node C by superscript *'C”, so 

M and K are subsets of M and K respectively. However, in actual computer 
implementation this element group is not necessary. We further denote M « 
and ]^= K^+ . 

Similarly, all vectors are. then partitioned accordingly into "A” and”B” 
parts, cf. 0 = V = and F * • The superscript ”T" 

denotes the tr^nspo^. The ^vector is sojpetimes r^efined by augmented 
matrices, £ = £ + £ where 0 * (£ £) and 0 (0 0 ^ . Similar defini- 
tions are jyged for V and F . Any nonzero terms in F obtairied in a compu- 

tation of F are neglected; they are assumed to be zero. 


As an example, consider a one dimensional mesh depicted in figure 1, it 
consists of 8 nodes and 7 elements. Then the set of nodes *’A” will be 1,2,3, 
4,5; the set of nodes will be 6,7,8; and the set of nodes "C” will be 5. 


ELEMENT 6 = 
NODE L 


© (D (D © I ® ® 0 

2 3 4 5 6 7 8 

I 

GROUP A (mAt) i GROUP B (At) 


Figure I. 
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Let M , K and F be the eth element mass, stiffness and force contributions 


to the global 

arrays 

respectively 

' ; then 
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*B 

If we let be the 1th component of the global assembled vector P 
1 ^ 
then ^ 

P =* ^Oj 0)0>0,P^ = 0,Pg >Py fPg) 

with these definitions ^ the mixed time partition Is given as follows. 
• Governing equation 
for j=0,m; 


MV + K*^e*^. + K®8*® = F . . 


(3.1) 


and 


for 


B *B B-*B C-'*A *B X 

M V ^. + K0 ^. + K 0 _ = F , . (3.2) 

~ ~n+j ~ ~n+j ~ ~n+j ~n+j 

^ / ^X <^.3C 

where 6 . - Is a suitable extrapolator (and/or Interpolator) of 9 (and/or 6 ) 

for x=A^aAd B. In actual computation, equation (3.2) Is ImpllciPly IncludeS^"*^ 
In equation (3.1), and for j=l,..*,m-l no quantities of A are being solved. A 
family of Integration partitions can then be deduced from equations (3.1) and 
(3.2) If M is assumed to be lumped. Some members which are of practical im- 

portances are shown in table 1. 


Table 1 


Designation 

Time Integration In 
Element Group A 

Time Integration In 
Element Group B 

Extrapolator/ Interpolator 

Node A 

Node B 

E-E 

explicit with At 

explicit with At 

0*^ 

'^n 

0*» 

mE-E 

explicit with mAt 

explicit with At 

0*^ 

~n 

0*® 

'^n 

mE-I 

Implicit with mAt 

explicit with At 

~n”hn 

0*® 

'^n 

E-I 

implicit with At 

explicit with At 

e*A 

~n+l 

0*® 

^n 

I-I 

implicit with At 

implicit with At 

0*^ 

'>^+1 

0*® 

~n+l 


For purposes of describing the computer Implementation and stability 
analysis, the modified generalized trapezoidal rule will be used to carry out 
the time temporary discretization of equations (3.1) and (3.2) though other 
implicit integration methods can also be used; they are: 


• Modified generalized trapezoidal rule 
f or j“l, . . . ,m; 



0^ + (l-a)jAt 
~n ~n 

for 1 < j < m define the set *’C” only. 


(3.3) 
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0® . . “ 0® . . 1 + (l-a)At V® . . T 
Ml+j '^+j “1 '^+j “1 

0^ + am At 


'^+m 


(3.4) 

(3.5) 


and 




0® 

Ml+j 


+ aAt 




(3.6) 


In the above equations, a Is a free parameter which governs the stability and 
accuracy of the method. We now Illustrate some useful partitions which have 
been depicted in table 1. 


Example 1 ; E-E partition 

In this case, m-1, for x«A and B. Equations (3.1) and (3.2) 

reduce to: 


MV ^, + K0 ^ = 


F 

^+1 


•^+i~ 


0 + (l-a)At V 
/N/n Ml 


(3.7) 

(3.8) 


and 


0 0 aAt V (3.9) 

~n+l ~n+l ~n+l 

Equations (3.7) to (3.9) represent the predictor-corrector explicit 
algorithms with equation (3.8) as the predictor and equation (3.9) as the 
corrector. 


Example 2 : mE-E partition 

In this case, m > 1, 
(3.6) reduce to: 



~n+j 


and 


0 ® 

~n+j 


0 ® 

-Sn+j 


Equations (3.1) to 


PREDICTOR PHASE: 


equation (3.3) (3.10) 

and 

equation (3.4) (3.11) 

GOVERNING EQUATIONS: 

equation (3.1) (3.12) 

and 

equation (3.2) (3.13) 
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CORRECTOR PHASE: 


equation (3.5) 


(3.14) 


and 


equation (3.6) 


(3.15) 


Example 3 mE-I partition ^ 

In this ca^y m > 1 ; In equation (3.1), i0 . . = 0 . for element group 
A only, portion which Is^’rAate^^ to K for j=0 and m. 

This Is done automatically If element group A Is defined to be the Implicit 
element group and element group^B Is defined to be the explicit element group. 
0* = 0 . for 1 < j < m, and 0 ■ ^ = 9 . , for 1 < j < m . Equations (3.1) to 

t^to) r^uce to: ~n j ~n j. 


PREDICTOR PHASE: 

equation (3.3) 

and 

equation (3.4) 
GOVERNING EQUATIONS: 
for j=0,m; 


(3.16) 


(3.17) 


MV . . + K^0*1, + . . = F . . 


(3.18) 


and 


for J=l,...,m-1; 


~ -:n+j 

CORRECTOR PHASE 


f K^e*® 


C-*A *B 
. + K^0 = F 

'^4*j Mi+j *^+3 


equation (3, 5) 
equation (3.6) 


(3.19) 

(3.20) 

(3.21) 


Example 4: E-I partition 

This is a special case of example 3. Equations (3.1) to (3.6) reduce to: 


MV + K^0 _ + K®0 =F 

'^H-l ^ ^ M [ l +1 Ml +1 

(3.22) 

0^=0 + (l-a)4t V 

(3.23) 

^+1 ^ ^+1 ^+1 

(3.24) 
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Equations (3.22) to (3.24) represent the Implicit-explicit algorithms devel- 
oped by Hughes and Liu (see e.g., [8,10]) in which equation (3.23) is the 

predictor and equation (3.24) is the corrector. 


Exas2^! I-I partition 

In this cose, m-I, s ^ ^ and 6^^, s 

reduce to the usual implicit formulation an<r it is: 


Equations (3.1) to (3.6) 


^+1 ^^+1 


-n+1 


-n+1 


e + (l-a)At V 
~n ~n 


-n+1 


6 . , + o4t V , , 

^^+1 Ml +1 


4 . STABILITY CRITERION 


(3.25) 

(3.26) 

(3.27) 


Our aim in this section is to deduce the stability characteristic of 
these mixed time partition algorithms. It suffices to restrict ourselves to 
the case in which F=0 and all mass matrices are lumped for purposes of sta- 
bility analysis. An energy balance technique (see [8J for a discussion) is 
employed to carry out the stability analysis. To simplify the subsequent 
writing, the following notations will be used. 


Fx , 1 = X . - X (4.1) 

'^-hn ^n 

<x > = [x + X )/2 (4.2) 

^n-hn '"^n+m ^n^ 

[x . 1 = X ...1 - X (4.3) 

and for j=0,l, . . . ,m-l 

<x ,.> = [x + X ,.1/2 (4.4) 

--n+j '^n+j +1 ~n+j ^ 


We have not made a complete satisfactory stability analysis of these mixed 
time partition procedures. However, if we assume: 


1. 

<V K <V , >/z K z > 1/m where 

*B *B m-2 *B 

s - Xn * j£l <Xn±J> 

<V*®>/V*^ > (1 +(m-l)a)^ 

~n ^ MU ~n ~ ~n 

(4.5) 

2. 

(4.6) 

3. 

<v*®, <V*®. >/V*^ V*S (1-a)^ for 

Mi+j -- --n 



J “1 > • • • >ni“l 

(4.7) 

and let 



4. 

T T 

*A *R *A *B *B *B 

S - V M V +V m"v >0 

n *^n Ml ^n '^n 

(4.8) 

5. 

- <v*^ ° 

n+m Mi+m Mi+m 

(4.9) 
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6 


(4.10) 


„B 


«*“>■' k" <V*?.> > 0 


the energy expression of these mixed time partition procedures can be shown to 
be: 

S ^ < S - 2mAt - 2At “zj; P®. . (4.11) 

_ n+m n n+m j =0 n+j 

R B C *R *B 

Here, K = K - K and the stability is governed by M and M provided 

a > 1/2 . Let 


and 


j3j - - l/2jAt 


Wj = ^ + (a - l/2)jAt 


(4.12) 


(4.13) 


the definitions of M and M for the five cases discussed in section 3 are; 


Example 1; E-E partition 




*B B 
and M = 


Example 2; mE-E partition 


M = Q1 


*B B 

^ =2i 


Example 3; mE-I partition 


*R C *B B 

M = + 0 and M = < 9 ^ 


'MU “^m 
Example 4; E-I partition 


M**^= ut + Q? 


:i ■ Si 
Example 5 : I-I partition 


*B B 
and M =5^ 


W? and M*^= W® 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


(4.18) 

*R 


J^ese mixed time partition procedures are stable if a > 1/2 and M 
and M are both positive definite. A summary of the results is as follows: 


Example 1 : E-E partition 

= 0® . < 2 
crit crit 

Example 2 : mE-E partition 

< 2 

crit crit 


and Q j ^ < 2 
crit 


(4.19) 


(4.20) 
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Example 3; mE-I partition 


«crlt ' 2 and 0 “^^^ < 2 


Example 4; E-I partition 




crit 

Example 5t I-I partition 


< 2 and 


unconditionally stable 

In equations (4.19) to (4.22), is defined to be JAt 
where denotes a typical eigenvalue of the eigenproblem 

M^0 + K*e = 0 


(4.21) 


(4.22) 


(4.23) 


(4.24) 


5. IMPLEMENTATION ASPECTS 


In this section, we generalize the mixed time integration methods de- 
scribed in section 3 to NUMEG element groups. Different time integration 

methods (implicit/explicit) with different time steps can be used in each ele- 
ment group. Let Atj,„^ and be the element group time step and element 

group time respectively for NEG = 1,..., NUMEG. There are NUMEL elements in 

each element group. We denote At as the minimum time step amount for all these 
element groups. In this formulation all element group time steps are required 
to be integer multiples of At and the time steps for adjacent groups are 
integer multiples of each other. Furthermore, no two implicit groups with 
different time steps can be adjacent to each other. In addition, for each 
implicit group that element group time step must be greater than those of the 
adjacent explicit groups. The main advantage of this m^ implicit - m 2 
explicit - m^ implicit -... etc. technique is to minimize the semi -bandwidth 
of complicated problems especially in the three-dimensional case. To illus- 
trate the idea, consider the one dimensional mesh shown in figure 2. It con- 
sists of NUMEG element groups and NUMNP nodes. In this case NUMEG is equal to 

4 and NUMNP is equal to 12. We assumed that node 1 is an essential boundary 

condition node and hence the number of equations, NEQ, is equal to 11. The 

essence of the present development can be deduced graphically by considering 


NODE 

GROUP 


I{6At) 


E(2A-t) 


I (4 At) 


E(At) 


I 23456709 


10 II 12 
4 


Figure 2. 
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the solution procedures of the matrix equations. The "active column equation 
solver” is the key to the success of this technique (see [8, 13] for a descrip- 
tion of this equation solver). The profile of the effective stiffness matrix 
K of this one dimensional mesh is shown in figure 3. We can observe from 
figure 3 the following: 

Group 1: implicit with At^= 6At , five words of storage (1-5), 3 elements and 
3 equations. 

Group 2: explicit with At 2 = 2At , two words of storage (6-7), 3 elements and 2 
equations . 

Group 3: implicit with At^= ^At , seven words of storage (8-14), 3 elements 

and 4 equations. 

Group 4: explicit with At^= At , two words of storage (15-16), 2 elements and 
2 equations. 


NODE NO 
2 

3 

4 

5 

6 

7 

8 

9 

10 
II 

12 * 


Figure 3. 



The equation systems of each element group are uncoupled and hence each group 
can be integrated at its owq^ group time step. For example, we assume the 
effective stiffness matrix K is formed and factorized once. In a time 
interval of 6At , group 1 will be integrated Implicitly once, group 2 will be 
integrated explicitly three times, group 3 will be integrated implicitly once 
and group 4 will be integrated explicitly six times. In order to handle the 
forward, reduction and backsubstitution and update procedures automatically, we 
required two arrays At „„ and J each has a dimension of NUMNP. At 

array contains the nodal ^ime steps of each node. Nodes associated with only 
one element group NEG are assigned a time step of At^^ , whereas those which 
are in common to other element groups are assigned to have the maximum time 
step from the adjacent groups. array contains the nodal time of each 

node. From these two arrays (At^^j^^ and the boundary condition 
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codes, another time step array 
equation systems can then be generate 


and an equation time array T of the 

equation systems can then be gene*^rated. Both At„„Q and ^MgQ have oimensions 
of NEQ. We require further a master time which is incremented by the 
smallest time step At . For this particular example the and Atj^g^ 

arrays are: 


At„^„„=r 6At , 6At , 6At , 6At , 2At , 2At , 4A t , 4At , 4At , 4At , At , At ) 
NODb ^ 


and 


AtNEQ=®( 6At , 6At , 6At , 2At , 2At , 4At , 4At , 4At , 4At , At , At) 

The T _ and arrays are incremented by time steps of At^^ g and 

At™ respectively. With these definitions, the generalized mixed time inte- 
graiSon is to proceed over the time interval . The procedures are as 

follows: 


1 . 


2 . 


3. 


where 


and 


4. 


5. 

6 . 


6a. 


Initialization 

Set Tjjgg,Tj^Qjjg and Tj^g^= 0 

Determine V 
~o 

V = m"^[F - K0 ) 

'k 

Form and factorize K 

* NUMEO *NEG 
K = I K 
~ NEG=1 ~ 


•NEG. if implicit 

/N/ ^ 


K*^EG= g diagonal matrix if explicit 

± * 

At : set effective force F equal to zero 

MM 

Loop on element groups NEG=^1, • • • ,NUMEG 
“ ^NEG+ ^*^NEG > ^a 

Loop on elements e=l, • • • ,NUMEL 

Define predictor values 0 


\0DE ^^NODE < 

Sode" -node'*' ^^"“^^^node ^ode 

^NODE '*' ^*^NODE ^ 


(5.1) 


(5.2) 


(5.3) 


(5.4) 


(5.5) 
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~NODE “ ~NODE '"■N0DE^~N0DE 

*e 

6b* Form element effective force f 

f*®= M®0® + aAt„„_F® if Implicit 

and 

f*®= aW(F®- K®0®} if explicit 

where 


(5.6) 


(5.7) 


(5.8) 


6c. 

6d. 

5a. 

7. 

7a. 


7b. 

8 . 

8a. 


8b. 

9. 


0 

W = diagonal matrix with along the diagonals 

Sum up effective force from element contributions 


End of element loop 
End of element group loop 

Solve for 9 , i*e., forward reduction and backsubstitution 


Loop on equation number N=1,...,NEQ 

^neq'*’ ^*^neq ^ \ 

Forward reduction and backsubstitution for equation N 
End of equation number loop 
Update V and 0 


Loop on N=l, • • . ,NUMNP 
^NODE'*' ^*^N0DE ^ 

y ^ (e - 0 


go to 8b 
7 


End of nodal number loop 

If Tjj > stop, otherwise go to 4 


(5.9) 


(5.10) 


(5.11) 


+ 'V* means "is replaced by" 


174 



6. NUMERICAL EXAMPLES 


The stability and accuracy of these mixed time partition procedures are 
confirmed by the following one dimensional heat conduction problem* The 
finite element mesh consists of a rod kept at a temperature of 0*0 for all 
time at the left end and insulated at the right end, and subjected to a con- 
stant initial teDq>erature of 0*1* The finite element model consists of (from 
left to right) 10 elements, each with a length = 10.0, 10 elements, each 
with a length of 100*0, and 10 elements, each with a length of 10*0* The 
thermal diffusivity, C^, is set to 200.0* A total of five computer runs are 
being made* They are (1) explicit, (2) explicit-implicit, (3) 10 explicit- 
explicit , ( 4) 10 explicit-implicit , and ( 5) implicit • All analyses are run 
with At ■ 0.25 which is the critical time step based on 10.0. For the m 
explicit-implicit cases, the implicit elements are the middla 10 elements 

100.0). The results obtained from these five analyses are compared to the 
analytical solution. They are depicted in figure 4* The temperature reported 
is at X = 100.0. The solution time ratios (normalized by the implicit time) 
for the above five cases are 0.892, 0.971, 0.578, 0.638 and 1.000 respec- 
tively. 




2 

In order to demonstrate the advantages of this mixed time technique C 
is raised to 2.0 x 10 for the implicit elements. A time step of At =* 0.25 
is used for the explicit-implicit and 10 explicit-implicit runs. This problem 
would therefore not be stable with a lOE-E partition. The results are pre- 
sented in figure 5. 

All the above calculations performed with a=0.5 and lumped mass matrices 
are assumed throughout. 
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7, CONCLUSION 


In this paper, we have developed a family of mixed time partition proce- 
dures for transient thermal analysis of structures. Both the stability cri- 
terion and the Implementaton aspects of these methods are described. Numer- 
ical corroboration of the stability and accuracy of these techniques is also 
presented. The implementation procedures of these new algorithms are straight- 
forward and are recommended for Inclusion in current thermal analysis computer 
programs. 
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